Thickness of the strangelet-crystal crust of a strange star 



OO 
O 

o 
< 

in 



p 

00 

o 

00 

o 



X 



Mark G. Alford and David A. Eby 
Physics Department, Washington University, St. Louis, MO 63130, USA 
(Dated: 5 August 2008) 

It has recently been pointed out that if the surface tension of quark matter is low enough, the 
surface of a strange star will be a crust consisting of a crystal of charged strangelets in a neutralizing 
background of electrons. This affects the behavior of the surface, and must be taken into account in 
efforts to observationally rule out strange stars. We calculate the thickness of this "mixed phase" 
crust, taking into account the effects of surface tension and Debye screening of electric charge. 
Our calculation uses a generic parametrization of the equation of state of quark matter. For a 
reasonable range of quark matter equations of state, and surface tension of order a few MeV/fm^, 
we find that the preferred crystal structure always involves spherical strangelets, not rods or slabs 
of quark matter. We find that for a star of radius 10 km and mass l.SM©, the strangelet-crystal 
crust can be from zero to hundreds of meters thick, the thickness being greater when the strange 
quark is heavier, and the surface tension is smaller. For smaller quark stars the crust will be even 
thicker. 

PACS numbers: 25.75.Nq, 26.60.+C, 97.60.Jd 



I. INTRODUCTION 

Quarks in their most familiar form are confined in pro- 
tons and neutrons that make up standard nuclear mat- 
ter. However, according to the "strange matter hypoth- 
esis" [U, [3j this form of matter might be metastable and 
the fully stable state would then be "strange matter", 
which contains roughly equal numbers of up, down, and 
strange quarks. Large (kilometer-sized) pieces of strange 
matter are "strange stars" (for a review see Ref. 0]); 
small nuggets of strange matter are "strangelets" Q • The 
strange matter hypothesis remains a fascinating but un- 
proven conjecture. In this paper we will assume that it 
is correct, and investigate the structure of the crust of a 
strange star. 

The traditional picture of the surface of a strange star 
is a sharp interface, of thickness 1 Fermi. Below the 
interface lies quark matter, the top layer of which is posi- 
tively charged. Above the interface is a cloud of electrons, 
sustained by an electric field which could also support a 
thin nuclear matter crust in suspension above the quark 
matter [1, [1] , as long as the strange star is not too hot Q ■ 

However, if the surface tension a of the interface be- 
tween quark matter and the vacuum is small enough, the 
surface will take on a much more complicated structure. 
If (T is less than a critical value CTcrit then large lumps of 
strange matter are unstable against fission into smaller 
pieces @, 0]. As a result, the simple surface described 
in the previous paragraph is unstable, and is replaced 
by a mixed phase involving nuggets of positively-charged 
strange matter in a neutralizing background of electrons. 
It is reasonable to guess that the ground state is a regu- 
lar lattice, leading to a crust with a crystalline structure. 
(Note that this crystal is completely different from the 
Larkin-Ovchinnikov-Fulde-Ferrell phase of quark matter 
[13) [Hi; where the quark matter density is uniform, but 
the pairing gap varies in space.) Jaikumar, Reddy, and 
Steiner conjecture that the strangelet crystal crust will 



actually be a multi-layer structure of mixed phases, anal- 
ogous to the "nuclear pasta" phases that occur in mod- 
els of the inner crust of a conventional neutron star [l^ . 
At the outer edge of the crust, we expect a dilute low- 
pressure lattice of small strangelets in a degenerate gas of 
electrons. As we descend in to the star, the pressure rises, 
and the structure is modified (becoming denser, and per- 
haps changing to rods, slabs, cavities, etc). At a critical 
pressure pcrit the mixed phase is no longer stable, and 
there is a transition to uniform neutral quark matter. As 
one burrows deeper into the star the pressure continues 
to rise, and there may be other phase transitions between 
different phases of quark matter [l^l, but those will not 
concern us here. Note that in this scenario the strange 
star depends on gravity for its existence. In the absence 
of gravity, it would undergo fission into strangelets. 

Ref. [8^, assuming zero surface tension and neglecting 
Debye screening, estimated that the mixed phase crust 
might be 40 - 100 m thick, with pcrit ~ 1000 MeV^. This 
is an interesting result because if a strange star has a suf- 
ficiently thick crystalline crust, it might be hard to dis- 
tinguish from the crust of a neutron star. Astrophysical 
properties that are sensitive to the crust include cooling 
behavior, neutrino and photon opacity during a super- 
nova, the photon emission spectrum, glitches, and fre- 
quencies of seismic vibrations which are observed after 
giant flares in magnetars. For further discussion and ref- 
erences see Sec. IVII This paper will make a more careful 
calculation of the properties of a strangelet crystal crust, 
including the effects of Debye screening [l3| and surface 
tension. 

We expect that the properties of the crust will emerge 
from a competition between various different contribu- 
tions to the energy. Charge separation is often favored 
by the internal energy of the phases involved, because a 
neutral phase is always a maximum of the free energy 
with respect to the electrostatic potential (see [H, IT5| : 
for a pedagogical discussion see [ia|). The domain struc- 
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ture is determined by surface tension (which favors large 
domains) and electric field energy (which favors small 
domains). Debye screening is important because it redis- 
tributes the electric charge, concentrating it in the outer 
part of the quark matter domains and the inner part of 
the surrounding vacuum, and thereby modifying the in- 
ternal energy and electrostatic energy contributions. 

To make an estimate of the thickness of the crust we 
need to calculate the equation of state of the mixed phase, 
i.e. the energy density e^np as a function of the pressure 
Pmp- The thickness of the crust for a star of mass M is 
then 

AR = RstJ^ - 2) — dp„,p , (1) 

in ?i = c = 1 units. This expression follows from the 
Tolman Oppenheimer Volkoff equation [3 ) assum- 
ing that AR <C i?star, and that everywhere in the crust 
the pressure is much smaller than both the local energy 
density and the average energy density of the whole star. 
These are very good approximations for the cases that 
we study. 

We obtain function of p^p by dividing the 

strangelet lattice into unit cells ( "Wigner-Seitz cells") 
and calculating the pressure at the edge of a cell as a 
function of its energy density. We study cells that are 
three-dimensional (a lattice of strangelets in a degenerate 
gas of electrons), two-dimensional (rods of strange mat- 
ter in a degenerate gas of electrons) and one-dimensional 
(slabs of strange matter interleaved with regions of de- 
generate electron gas). Our approach is similar to that 
used in studying mixed phases of quark matter and nu- 
clear matter in the interior of neutron stars [Toj . 

We build on the formalism for a generic quark matter 
equation of state and infinitely-large Wigner-Seitz cells 
that was developed in [sl, . The main assumptions that 
we make are: 

1) Within each Wigner-Seitz cell we use a Thomas-Fermi 
approach, solving the Poisson equation to obtain the 
charge distribution, energy density, and pressure. This 
is incorrect for very small strangelets, where the energy 
level structure of the quarks becomes important 20, 2l|; 
such corrections may be relevant for the very low pres- 
sure (outer crust) part of our results (see Sec. IVip. 

2) We assume our D-dimensional Wigner-Seitz cells to be 
Z?-spheres. In reality the cells will be unit cells of some 
regular lattice (cubic, hexagonal close packed, etc). How- 
ever, as long as the cell is much bigger than the strangelet 
inside it, we expect this approximation to be reasonably 
accurate. We will only report results for cases where 
i?ccii > 2i? (i? being the radius of the quark matter in the 
center of the cell, which we expect will have a rotation- 
ally symmetric shape because of the surface tension) . In 
some cases this assumption is violated, and we will then 
only be able to obtain a lower limit on the crust thickness 
(see Sec. lVBl]) . 

3) We treat the interface between quark matter and the 
vacuum as a sharp interface, with no charge localized on 



it, which is characterized by a surface tension. We neglect 
any surface charge that might arise from the reduction 
of the density of states of strange quarks at the surface 
[2^ m, 113, . We also neglect the curvature energy of 
a quark matter surface [2^, [l^l , so we do not allow for 
"Swiss-cheese" mixed phases, in which the outer part of 
the Wigner-Seitz cell is filled with quark matter, with a 
cylindrical or spherical cavity in the center, for which the 
curvature energy is crucial. Note that these phases would 
be expected to occur at higher pressure than the ones we 
study, so including them is likely to make the crust even 
thicker than we predict. 

4) We assume that the chemical potential for negative 
electric charge fie is much less than the chemical potential 
for quark number /i. This allows us to expand the quark 
matter equation of state in powers of /ie, and means that 
within the quark matter we can ignore the contribution of 
electrons to the charge and pressure. This is a very good 
approximation for small strange quark mass, which cor- 
responds to small Uq in our parameterization. For the 
largest value of Uq that we study, /ig in neutral quark 
matter is close to 100 MeV, and the assumption is still 
reasonable. 

5) We assume that only electrons are present, with no 
muons. This is valid as long as /Xe is less than the muon 
mass m^, which is true for all the cases that we study. 

6) We assume that /ig is always much greater than the 
electron mass. Thus in the degenerate electron gas, we 
can take the electrons to be massless, which simplifies 
the Thomas-Fermi calculation of their charge distribu- 
tion. Since fj,e drops monotonically from the center of 
the cell to its edge, this condition will only be violated 
for very large cells (very low pressures). 

8) We always work at zero temperature. The tempera- 
ture of the surface of a compact star, even during a flare 
is expected to be less than 100 keV, so we expect 
this to be a reasonable approximation. 

II. CHARACTERIZATION OF QUARK 
MATTER 

A. Generic parametrization 

We use the generic parametrization of the quark mat- 
ter equation of state suggested in Ref. 

PqmIa*- Me) = Po(m) - + \Xq{p)pI + • ■ ■ (2) 

which expresses the pressure as a function of the chemical 
potential for quark number (/i) and for negative electric 
charge (/Xe), expanded to second order in /Xe. In addition, 
we assume that there is a surface tension a associated 
with the interface between quark matter and vacuum. In 
this paper we do not include curvature energy. 

This parametrization is model- independent. Any spe- 
cific model of quark matter can be represented by appro- 
priate choices of a, po, the charge density rig, and charge 
susceptibflity Xn- 
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The quark density n and the electric charge density 
'^'qm (^'^ units of the positron charge) are 



dp, 



QM 



9qm 



dp, 



QM 



So in uniform neutral quark matter the electron chem- 
ical potential is ^n<="t'ai = "q/Xq- Eq. ([2]) is a generic 
parametrization if ^nout'^i <^ which is typically the 
case in three-flavor quark matter. 

The bag constant enters in poilJ-), and we will fix it by 
requiring that the first-order transition between neutral 
quark matter and the vacuum occur at quark chemical 
potential ^crit, i-e. p(/iciit, Me''"*'^''') — 0- Because we are 
assuming that the strange matter hypothesis is valid, we 
must have /icrit ^ 310 MeV, since at /i w 310 MeV there 
is a transition from vacuum to neutral nuclear matter. 
In this paper we will typically use /icrit = 300 MeV. The 
value of fj, inside our quark matter lumps will always be 
very close to /icrit, so we can evaluate Uq and Xq at /icrit, 
and not be concerned about their //-dependence. 

We will restrict ourselves to values of the surface ten- 
sion that are below the critical value 191 



dcrit = 0.1325- 



XQ 



0.1325- 



3/2 ' 



where a = 1/137 and Au is the Debye length 

1 



A 



D 



(4) 



(5) 



For typical models of quark matter, (Tcrit is of order 1 to 
10 MeV/fm^ (see Table H]). If the surface tension is larger 
than acrit then the energetically favored structure for the 
crust will not be a strangelet crystal but the simple sharp 
surface that has has been assumed in the past [g, 0] ■ 



B. Specific equations of state 

When we show numerical results we will need to vary 
Tig and Xq over a range of physically reasonable val- 
ues. To give a rough idea of what values are appropriate, 
we consider the example of non-interacting three-flavor 
quark matter, for which Uq and Xq become functions of 
/i and the strange quark mass rus , while po is in addition 
a function of the bag constant B. Expanding to lowest 
non-trivial order in m^. 



nQ{fi,ms) 
XQ(/i,m,) 



27r2 



B 



(6) 



We emphasize that these expressions are simply meant 
to give a rough idea of reasonable physical values for Uq 



and Xq ■ Our treatment does not depend on an expansion 
in powers of nrig. To tune the transition between neutral 
quark matter and the vacuum so it occurs at /i = /icrit 
(see previous subsection), we set B so that po (Merit) = 

^ riQ (Merit )/Xq (/icrit)- 

In the regions between lumps of strange matter, we will 
assume that there is a degenerate electron gas, which we 
treat in the Thomas-Fermi approximation. As long as 
/ie is much greater than the electron mass, we can treat 
the electrons as massless particles, whose pressure and 
charge density (in units of e) is 



Pe~ ifJ-e) 



127r2 



ge-(Me) = -ITT • 



(7) 



III. ANALYSIS OF A WIGNER-SEITZ CELL 

We will study one, two, and three dimensional Wigner- 
Seitz cells. In the center there is a slab, rod, or sphere of 
quark matter, with radius R. The cell itself has radius 
^ccii- We want to calculate the equation of state of a 
mixed phase made of Wigner-Seitz cells, so we solve for 
the charge density, energy density, and pressure through- 
out the cell, using the Thomas-Fermi approximation for 
the contributions of quarks and electrons. This corre- 
sponds to solving the the Poisson equation, which reads 
(in Heaviside-Lorentz units with h = c = 1) 



(8) 



where q{r) is the electric charge density in units of the 
positron charge e, and fie is the electrostatic potential 
divided by e. The equation is not trivial to solve because 
the charge density is itself a function of fie (see (jS])). The 
boundary conditions are that there is no electric flcld in 
the center of the cell (no (5- function charge there) , and no 
electric field at the edge of the cell (the cell is electrically 
neutral). 



^(0) = 0, ^(i?ccn) = 
ar ar 



(9) 



We also need a matching condition at the edge of the 
quark matter, i.e. at r = i?. As discussed in Sec. HI we 
assume that no charge localized on the surface, so we 
require continuity of the potential and electric field at 

fieiR-i^) = fieiR-S) , ^{R^) - ^(R-S) . (10) 

ar ar 

In two or three dimensional cells, the value of jj, inside 
the strange matter will be slightly different from /icrit 
because the surface tension compresses the droplet. To 
determine the value of /i, we require the pressure discon- 
tinuity across the surface of the strangelet to be balanced 
by the surface tension: 



PQM(M,Me(i?)) -Pe-(Me(^)) 



R 



(11) 
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Once these equations are solved, we can obtain the 
relevant properties of the cell. The total energy of a D- 
dimensional Wigner-Seitz cell is 



as 



E 



+ 



R 



(12) 

where fie is a function of r, and floir) is the surface areas 
of a (D — l)-sphere, i.e. 



ill 



(r) = 27rr , n^ir) = Airr"^ 



(13) 



The —^iJeQ terms in (fT^ come from combining —fieQ 
(from the relationship between energy density and pres- 
sure) with the electric field energy density +^fieq- The 
external pressure of the cell is simply the pressure of the 
electrons at the edge of the cell, 



PcKt =Pe-(Me(-Rcoll)) 



The total number of quarks is 



N 



flDir)drn{fi,fie) 



(14) 



(15) 



The volume of the cell is ^ = 2Rce\i, nRceii , or 
(4/3)7r_Rccii'^ for D = 1,2,3 respectively, because we are 
assuming rotationally symmetric cells. 

By varying R and i?cGii we generate a two-parameter 
family of strangelets. However, there is really only a 
single-parameter family of physical configurations, pa- 
rameterized by the external pressure Pcxt- On each line of 
constant pext in the (i?, i?ccii) parameter space, we must 
minimize the enthalpy per quark. 



h = 



E + Pe^tV 

N 



(16) 



to find the favored configuration. We are at zero temper- 
ature so h is also the Gibbs free energy per quark. This 
is done separately for £> = 1, 2, 3 cells, and the structure 
with the lowest h is the favored one. 

We now have a well-defined way to obtain the equation 
of state of the mixed phase of quark matter, namely the 
energy density £mp = E/V as a function of the pressure 
Pmp = Pcxt- This, via ([1]), determines the thickness of the 
strangelet crust. 



IV. SOLUTIONS FOR THE WIGNER-SEITZ 
CELL 

A. Quark matter 

Inside the quark matter, we can solve the Poisson equa- 
tion ([5]) analytically. We can rewrite it using ^ and ([3]) 



V^^e(r) = -47ra(nQ - XQMe(»'))- 



(17) 



In _D = 1, 2, or 3, the solutions obeying the first boundary 
condition in ^ are 



^J■e,lD[r) = hAcosh( — ) , 

XQ Ad 

Me,2D r = ^+AJo — , 
XQ Az3 

^J'e.3D{r) = h- smh( — , 

XQ >^Dr \d 



(18) 



The function Jg is the zeroth order Bessel function of 
the first kind, and it is a function of the square of its 
argument, so the result is always real. The integration 
constant A will be determined by matching to the solu- 
tion outside the strange matter. 



B. Electron gas 

In the degenerate electron gas region outside the 
strange matter, from ([7]) and ([5]) the Poisson equation 
becomes 



(19) 



There are three ways in which we were able to solve this 
equation. In D = 1, there is an exact analytic solution, 
which we present below. In any number of dimensions 
there is an approximate analytic solution, obtained by 
perturbing in powers of a, which works as long as the 
cell is not too large. Finally, one can use brute-force nu- 
merical methods to solve the differential equations with 
the appropriate boundary conditions. We used all three 
methods, checking their agreement with each other in 
situations where more than one of them was applicable. 
In our numerical results we give the values obtained by 
numerical solution of the Poisson equation. 



1. Analytic solution for slabs 

In one dimension, the Poisson equation is 
d^/ie 4a , 



(20) 



By a change of variable to (p = i-^2a/37r/Lte, this becomes 

-2^3 ^ (21) 



which belongs to a class of differential equations whose 
solutions are Jacobi elliptic functions sn(r|m) where m is 
the "parameter" [2^ . (Some authors write this as sn(r, k) 
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where m = k"^ and k is the elhptic modulus.) The Jacobi 
elliptic function obeys 



(Psn{r\m) 



(1 + m) sn(r|m) + 2msn{r\m) 



(22) 



which reduces to (PT|) for m = —1, so the closed form 
solution is 



ip{r) = — sn(^ — + iK{-l) , 



(23) 



where X is an integration constant, to be fixed by match- 
ing to the quark matter solution at r = i?. The argu- 
ment is shifted by iK{—l) to ensure that Lp'{Rcc\\) — 
(the boundary condition ([5]) at r = i?con)- K{m) is the 
complete elliptic integral of the first kind, so K{—1) k, 
1.3110288. Then ip is purely imaginary, and undoing the 
change of variables, the solution for /ig is 



The perturbative solution works when screening is a 
small correction to the unscreened (zeroth-order) electro- 
static potential. It is most reliable in three dimensions, 
where the zeroth-order electrostatic potential becomes 
small at large r. We used it to check our numerical re- 
sults. 



V. NUMERICAL RESULTS 

To get a good estimate of how thick strange star crusts 
might be, we vary /^crit, "-g, and Xq in the quark matter 
equation of state ([6]), and the surface tension cr, over a 
physically reasonable range, calculating the crust thick- 
ness in each case. The results of our calculations are 
displayed in Tables HI and HIl Before discussing them, we 
describe how they were obtained. 



t^e,iD[r) = -«Y^Y^°V X + 

(24) 

For one-dimensional slabs of quark matter we now have 
a complete analytic solution, combining ([18]) at r < i? 
with (I21D at r > i?, with X and A fixed by 



2. Perturbative solution for cylinders and spheres 



Another approach to solving ()19|) . which works in any 
number of dimensions, is to expand in powers of the elec- 
tromagnetic coupling strength a. We write 



/Xe(r) = /xo(r) + Q;/ii(r) + iJL2{r) + 



(25) 



Substituting in to and identifying powers of a, we 
find that to order a. 



V2^o(r) = 



vVM--/io(0' 

37r 



(26) 



Solving these equations in D = 2, 3, wc find 



Me,2D(r) = B + C\og{r/\D) + — [2B^ - &B'C 

+ 9BC^ - 6C^) + 3{2B^ - ABC + 3(7^) log(r/Ai3) 

+ 6{B - C)C^ logir/Xof + 2C^ log(r/AD)' 

l^eMr) =B-- + ^ \B\^i9C + Br) 
- 6C\C^Br)\og{r/XD) ■ 

(27) 

In each equation, the first two terms on the right-hand 
side are the vacuum solution /j,o, and the remainder are 
the first-order correction. The integration constants B 
and C, along with A from ([TS]). are determined by the 
boundary condition ^ at r = Rcdi and the matching 
condition ^TU\i . 



A. Geometry of the mixed phase 

For a given quark matter equation of state, we need to 
find the maximum pressure up to which a stable mixed 
phase exists. At each lower value of the pressure we must 
establish the geometric configuration of the mixed phase. 
We can then calculate the energy density as a function 
of the pressure, and obtain the crust thickness using ([T]). 

We follow the procedure described at the end of 
Sec, mil varying the radius of the cell -Rcoii and the radius 
R of the quark matter region at its center, to find the cell 
configuration with the lowest enthalpy per quark at each 
value of the pressure. In Fig. [T] we show some examples 
of the search for the favored cell configuration at a given 
pressure. We plot the "excess enthalpy per quark" 



Ah — h — hr. 



(28) 



as a function of R, for each of the three geometries, and 
for various different pressures. Here h is the enthalpy 
per quark (|16p in a given cell, and and ft-oo is its value 
in uniform neutral quark matter of the same pressure. 
Configurations with negative Ah therefore correspond to 
stable mixed phases at the given pressure. 

The results in Fig. [T] are for quark matter with ficnt = 
300 MeV, Uq = 0.0791 (corresponding to = 200 MeV 
in ([H])), and Ad — 4.82 fm (again, a value appropriate 
to free quark matter with fi = 300 MeV). The surface 
tension is 0.3 MeV/fm^. The first panel of Fig. [T] shows 
Ah{R) for 3D cells. The upper (red) curve is Ah{R) at 
the critical pressure PcHt = 1370 MeV^, which is defined 
by the presence of a minimum with Ah — (at R w 
10 fm in this case). We also show Ah{R) at a lower 
pressure, p — 10 MeV*; now there is a clearly favored 
mixed phase, with strange droplets of radius i? w 4 fm. If 
we push the pressure down to zero then the cell size goes 
to infinity, and the minimum in the Ah{R) curve moves 
further down to around —0.75 MeV. In the second panel 
we show Ah{R) for 2D cells at the same two pressures. It 
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Enthalpy per quark, Droplets (3D) Enthalpy per quark, Rods (2D) Enthalpy per quark. Slabs (ID) 




droplet radius (fm) rod radius (fm) slab radius (fm) 



FIG. 1: Search for a stable mixed phase at different pressures. We show the excess enthalpy per quark A/i (|28|) as a function 
of the radius R of the quark matter region in the center of the cell. At each value of R, the cell radius has been chosen to 
give the desired pressure. The favored configuration is the one with the smallest Ah, where Ah — for uniform neutral quark 
matter of the same pressure. For 3D cells, we find stable mixed phases (minima with negative Ah) for pressures below a critical 
value, pcrit = 1370 MeV*. 2D cells always had a higher Ah; ID cells were never stable. These plots are for a quark matter EoS 
with ^crit = 300 MeV, Xd = 4.82 fm, Uq — 0.0791 fm""^, corresponding to free quarks with a moderately heavy strange quark 
{ms = 200 MeV in dSJ). The surface tension was a = 0.3 MeVfm"^. 



is clear that the 2D structure has a lower critical pressure, 
and at these two pressures it is energetically unfavored 
relative to the 3D structure. In the third panel we show 
Ah{R) for ID cells. These appear to be even less favored. 
At p = 1370 MeV* the Ah{R) curve is already almost 
at its zero-pressure limit, which is never negative and 
therefore allows no mixed phase. We had to show Ah{R) 
for a higher pressure, p = 10^ MeV*, to see any change 
in the curve. We conclude that for this quark matter 
equation of state and surface tension, and at the pressures 
studied in Fig.[Tl the only mixed phase that occurs is the 
3D (droplet) one. 

We note the following features of the favored configu- 
ration of the Wigner-Seitz cell: 

• Increasing the pressure disfavors mixed phases: the 
Ah{R) curve rises and minima are smoothed out. We 
hypothesize that this is because the pressure is deter- 
mined by the value of iie at the edge of the cell (fH)) : 
if /Xe(i?ceii) is increased then, because iJ-eir) is mono- 
tonic, /Ze in the quark matter is also larger (closer to 
^neutral But thls decrcascs the energy benefit of mak- 
ing a mixed phase, which arises from the departure 
from neutrality. 

• As the dimensionality of the mixing geometry de- 
creases from 3 to 1, mixed phases become less favored 
(at least in this range of pressures). We hypothesize 
that this is because in lower dimensional structures, 
a smaller proportion of the quark matter is near the 
surface, where fie is different from ^neutral ^ 

Since it is only the minima of Ah(R) that are phys- 
ically important, we focus on them, and in Fig. [5] we 
plot the value of Ah at the minimum as a function of 
pressure. The vertical dashed lines mark the pressures 
p = 10 MeV* and p = 1370 MeV* used in Fig. [T] We see 
that, for the values of the quark matter parameters stud- 
ied in these figures, only droplet (3D) mixed phases will 
occur: slabs are never stable, and rods are never favored 



Excess enthalpy per quark 

H=300 MeV, Xd=4.82 fm, nQ=0.0791 fm-\ 0=0.3 MeV fm-2 




'■^ 10 100 1000 

pressure (MeV^) 



FIG. 2; Minimum value of Ah as a function of pressure, 
for quark matter with the same characteristics as in Fig. [T] 
Dashed lines mark the pressures (10 MeV* and 1370 MeV) 
that were used in Fig. [T] We see that for this type of quark 
matter the 3D (droplet) structure is always energetically fa- 
vored over rods. The slab structure is never stable. 



over droplets. 



B. Crust thickness calculation 

In Table U we show the results obtained by repeating 
the calculations described above for seven different quark 
matter equations of state and four values of the surface 
tension. In each case, once we have established the en- 
ergetically favored configuration of the mixed phase at 
each pressure, we obtain the crust thickness using ([1]). 
In order to present numbers whose physical interpreta- 
tion is clear, we assume that the quark star has radius 
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10 km and mass 1.5 M©, so we can give an explicit crust 
tliickness in meters. From the prefactor in Eq. ^ it is 
easy to find the multiplicative factors that convert our 
values into thicknesses of crusts on stars with different 
masses and radii (Table [TH) . 

The first two columns of Table H] specify the quark 
matter equation of state by giving the value of Uq 
and the value oi Xd (which fixes Xq via ([5])). We fix 
Merit — 300 MeV. The third column gives the maxi- 
mum surface tension dcrit for which a crust of droplets of 
strange matter could occur, the fourth column gives an 
estimated upper limit on the thickness, and the last four 
columns give the thickness of the crust for various values 
of the surface tension a. 



1. Range of parameters studied 

Our assumption that the strange matter hypothesis 
is valid requires that ^crit must be less than the quark 
chemical potential of nuclear matter, about 310 MeV, so 
we fix /^crit = 300 MeV. The value of /i inside our strange 
matter lumps will always be within a few MeV of /^crit, 
because in order to get any crust at all the surface tension 
cannot be large enough to cause significant compression. 

Typical values of Xq will be around 0.2//^^;^. ([5]), cor- 
responding to Xd K 4.82 fm. Ref. found that in the 
2SC phase Xq is smaller by a factor of 2. In Table |T] 
we explore this range, using three values, Xq = 0.2/i^j.;^, 
Xq = O.l/i^j.jt, and Xq = 0.05/1^^.;^, corresponding to 
Xd ^ 4.82 fm, Xd ^ 6.82 fm, and and An = 9.65 fm 
(via ©). 

Typical values of Uq will be around 0.05/icrit"^g (O, 
and a reasonable range would correspond to varying rris 
over its physically plausible range, from about 100 to 300 
MeV. (To have strange matter in the star, must be 
less than ^lcnt■) In Table [J we use rig = 0.0445, 0.0791, 
and 0.124 fm^'^, which would correspond to tUs = 150, 
200, and 250 MeV in dH). For A^, = 9.65 fm we only 
show results for Uq = 0.00445 fm^"^. We do not show re- 
sults for Uq ~ 0.0791 fm^^ or Hq ~ 0.124 fm~'^, because 
in those cases the value of ^^''"''■^i would be 133 MeV 
and 208 MeV respectively, which violates our assump- 
tion that ^-iid is also above the muon mass 
nifj, = 105.66 MeV, so we would have to take into ac- 
count muons as well as electrons. 

The value of the maximum surface tension Ccrit for 
which a crust of droplets of strange matter could occur 
(third column of Table [l| follows from It is the maxi- 
mum surface tension at which an isolated (zero-pressure) 
droplet would have lower enthalpy per quark than neu- 
tral quark matter at zero pressure, i.e. at the onset phase 
transition ai — /icrit- The last four columns of Table U 
give our results for the thickness of the crust at a range 
of values of the surface tension a. For values of a above 
CTcrit there is no crust. The values of a that we use are 
physically reasonable, given that rough estimates of sur- 



Droplet and cell size in crust 

|I=300 MeV, A,o=4.82 fm, nQ=0.0791 fm-\ a=0.3 MeV fm-2 
1000 b j 




-I20 -100 -80 -60 -40 -20 
position relative to star surface (m) 

FIG. 3: Size of the quark matter droplets and the Wigner- 
Seitz cell (in fm) as a function of position in the the 120m- 
thick crust that we predict for quark matter with the stated 
parameter values, which correspond to free quark matter with 
a moderately heavy strange quark {nis — 200 MeV in Q). 



face tension from the bag model are in the range 4 to 
10 MeV/fm^ 



2. Crust thickness results 

In Table [H we present the results of our calculation 
of the crust thickness as a function of the quark matter 
parameters. Clearly the thickness is very sensitive to the 
values of these parameters. 

In Table HI] we give the correction factors for some typ- 
ical masses and radii that are expected for quark stars 
(see, for example, Ref. Q, Fig. 28). We see that smaller 
and lighter quark stars have thicker crusts. In fact, for 
a star of radius 4 km the crust thickness could eas- 
ily be of order 1000 m, at which point the assumption 
Ai?/i?star used in Eq. ^ begins to become questionable. 
If the crust is thick enough then one must keep higher 
powers of Ai?/i?star in calculating it, and it may even be 
necessary to solve the full Tolman Oppenheimer Volkoff 
equation [13, El- 

Some of the crust thicknesses in Table H] are given as 
lower limits. This is because for large A^i and rig, and 
low surface tension, mixed phases are highly favored, 
and persist up to extremely high pressures (of order 
10^ MeV^). To achieve such a high pressure requires 
a very "cramped" cell geometry, in which i?ccn is only 
a little larger than the quark droplet radius R. In such 
a geometry we can no longer trust our approximation 
of treating the cell as spherical, rather than a unit cell 
of a crystal lattice (e.g. cubic). We therefore use an al- 
ternative upper limit Pmax on the pressure: Pmax is the 
pressure below which the favored cell configuration al- 
ways has i?ccii > 2i?. We can then calculate a lower limit 
on the crust thickness, by integrating ([1]) up to Pmax- 
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0.0791 


1.69 
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67 
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4.82 


0.124 


4.12 
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220 


160 


39 




6.82 


0.0445 


1.51 


72 


40 
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6.82 


0.0791 


4.8 


230 


> 170 


120 


45 




6.82 


0.124 


11.6 


550 


>460 


>390 


280 


39 


9.65 


0.0445 


4.27 


140 


110 


75 


22 





TABLE I: Crust thickness AR (in meters) for a strange star of radius 10 km and mass 1.5 Mq; for more general stars see 
Table |TI] We calculate the crust for seven different quark matter equations of state and four values of the surface tension. The 
first two columns, Ad and Uq, specify the quark matter equation of state ([2]) (via ([5])). The third column gives the maximum 
surface tension for which a strangelet crust will occur Q. The fourth column gives a rough upper limit on AR using an 
estimate from Ref. 0. The last four columns give our results for the crust thickness in meters for different values of the surface 
tension a (given in MeVfm"'^) of the interface between quark matter and vacuum. The crust thickness is very sensitive to the 
equation of state and the surface tension. It ranges from zero up to several hundred meters. 



R (km) M (Mq) Ai?/(A7? in TableHl) 

4 0.05 8.3 

6 0.2 4.4 

8 0.5 2.8 

10 1.0 2.0 

TABLE IL Correction factors to be applied to crust thick- 
ness values in Table |IJ for quark stars of various radii and 
masses. These follow from the factor multiplying the integral 
in Eq. ((TJ. Note that smaller stars have thicker crusts. 



izing background of electrons. We never found any pres- 
sure for any quark matter parameters where 2D (rod) 
or ID (slab) structures were energetically preferred. It 
would be interesting to see whether this remains true 
when the cell is allowed to be a different shape from the 
strangelet (e.g. square or cubic cells). 

3. Internal structure of the crust 



In the limit of low surface tension, our results are com- 
patible with the upper limit Ai?niax (fourth column in Ta- 
bic |T| which follows from taking the estimate obtained by 
ignoring surface tension and Debye screening in Ref. [8| 
(their Eq. (13)) and applying the TOV correction factor 
(1 — 2GM / Rstax)- We find that even a relatively mod- 
erate surface tension, around 1 to 10 MeVfm"^, reduces 
the thickness of the crust or even eliminates it completely. 
The mechanism is clear: larger values of the surface ten- 
sion disfavor the mixed phase by increasing surface en- 
ergy costs, leading to a lower pcrit and thinner crusts. 
However, it remains possible that a quark star could have 
a crust hundreds of meters thick. 

We find that the crust is thickest for large values of 
Tig and Ad, as one would expect from the estimate 
Ai?inax Uq / Xq oc "■q^A|, 3 • This is consistent with the 
observation from Fig.[l]that we can measure how favored 
a mixed phase is by the depth A/imin of the minimum in 
Ah{R). Since the Ah{R) curve moves up as the pressure 
increases (see Sec. IV A[) . we can guess that the deeper 
the minimum at zero pressure, the thicker the crust. 
From Ref. ^ (Fig. 2 and Eq. (26)), we find that at a 
given (7, the depth of the minimum is Ahmin oc ~'^q/Xq, 
i.e. A/imin oc — rtgA|). Hence the mixed phase is more 
favored, and the crust is thicker, for larger values of ng 
and Xo- 

All the crusts in Table U consist entirely of 3D struc- 
tures, i.e. spherical droplets of quark matter in a neutral- 



In Fig. [3] we select one of our quark matter equations 
of state, and show how the properties of the strangelet 
crystal lattice vary with position in the crust. The hor- 
izontal axis is Ar = r — i?star, so more negative values 
correspond to deeper parts of the crust. The plot should 
end at Ar = 0, but we were not able to push our nu- 
merical calculations to that value, so the curves end at 
Ar « -10 m. 

As one approaches the surface, the size i?,ccii of the 
Wigner-Seitz cell grows very large (note that Fig. [3] uses 
a logarithmic scale for the vertical axis). This means that 
the strangelet density becomes very low. We expect that 
as Ar 0, i?ceU will diverge, since the pressure must go 
to zero, so fie at the edge of a cell must go to zero, so the 
cell size must become infinite. In this limit the droplets in 
the crust efi^ectively become isolated strangelets, and we 
expect their size to settle down to that of the most stable 
isolated strangelet for this form of quark matter. We can 
predict this value from Ref. Q, eqn (24): by minimizing 
the free energy per quark Ag/n (which is equivalent to 
Ah in this paper) for the values of ficnt, ng, Xq, and a 
used in Fig. O we find that the most stable strangelet has 
a radius of 3.0 fm, which is exactly the asymptotic value 
emerging in Fig. [3) For such small strangelets we expect 
that our Thomas-Fermi approach is not accurate, and 
shell model corrections will become important: including 
such corrections is a topic for future research. 

As one goes deeper into the star (Ar becoming more 
negative), the pressure rises, so the cell size decreases, 
and the droplet size slowly increases, until we reach the 
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critical pressure at which uniform neutral quark matter 
becomes favored over the strangelet lattice. This is a 
first-order phase transition, as is clear from Fig. [U so 
the curves in Fig. [3] end suddenly, without any singular 
behavior. In this paper we do not take in to account the 
possibility of a metastable lattice that might persist to 
higher pressures. 

Fig.[3]shows that the strangelet crystal crust of a quark 
star tends to be fairly dilute: over most of the crust the 
quark matter droplet size is small, of order the Debye 
length in quark matter A^, while the cell size is larger, 
by a factor of 10 or more. 



VI. DISCUSSION 

The calculations described in this paper give us a more 
precise picture of the strangelet-crystal crust of a quark 
star. The results presented in Table |T] show that the 
thickness of the strangelet-crystal crust of a strange star 
is very sensitive to the surface tension <j of the inter- 
face between quark matter and the vacuum, and to the 
quark matter parameters tlq and Xq @, which deter- 
mine the response of the quark matter to deviations of 
the electrostatic potential from its neutrality value. Our 
results are compatible with those of Ref. Q where an up- 
per limit on the crust thickness was obtained by ignoring 
surface tension and Debye screening. The crust is thick- 
est for large rtg and small Xq (large Ad). As discussed 
in Sec. IVB 21 we find that values of surface tension in 
the physically expected range, around 1 to 10 MeVfm^^, 
reduce the thickness of the crust and may even eliminate 
it completely, but it remains possible that a quark star of 
radius 10 km could have a crust several hundred meters 
thick. From Table |TT] we see that for a smaller star the 
crust could be even thicker. 

The geometry of the mixed phase in our crusts, on 
the other hand, shows no variation at all. It is always 
three-dimensional, containing spherical droplets of quark 
matter. We never find any case where a two-dimensional 
(rod) or one-dimensional (slab) geometry is favored. 

Our calculations and results suggest two directions 
for future work: firstly, one could study phenomenologi- 
cal consequences of our understanding of the strangelet- 
crystal crust of a quark star. Secondly, one could improve 
on our treatment of the strangelet crystal, by relaxing 
some of the assumptions listed in Sec. ID 

The most obvious phenomenological task is to re- 
visit computations of the frequencies of seismic vibra- 
tions which are observed after giant flares in magnetars 
[s^ , [S^ l . Ref. [1^ found that the strangelet crystal crust 
did not have the right spectrum of toroidal shear modes 
to account for current observations: it would be interest- 
ing to see whether taking in to account the surface ten- 
sion and Debye screening affects that conclusion. Other 
aspects of the phenomenology of the crust could also be 
studied, for example (a) the thermal response of the crust 
to accretion [35i] ; (b) the role of the crust in the trapping 



of neutrinos and photons just after a type-II supernova 
f36|; (c) the spectrum of photons radiated from the sur- 
face of a quark star [s^, [H, [3§|; (d) the contribution 
of the crust to the moment of inertia and glitches ^4^; 
(e) the damping of r-modes in by shear viscosity in the 
crust [U, (for quark stars, the contribution from the 
interior has been calculated [H, IS]); (f) the thermal 
relaxation time of the crust and its response to the post- 
supernova "cooling wave" [i^. The thermal relaxation 
time of the crust depends on the thermal conductivity, 
for which we can make a very rough estimate using ap- 
pendix A of Ref. [i^. We find values of order a few 
hundred MeV^ at T 0.1 MeV, which is comparable to 
the range 10^^ ergcm~^s"^K~^ for low-density nuclear 
matter (Ref. [45] , Fig. 4) . We defer a more accurate cal- 
culation to future work. 

To improve on our treatment, the most pressing is- 
sues are to use a realistic shape for the Wigner-Seitz 
cells (which should be unit cells of some regular lattice, 
rather than spheres), to include shell- model corrections 
for the smallest strangelets, and to allow for "Swiss- 
cheese" phases where most of the unit cell consists of 
quark matter, with a hole at the center containing elec- 
trons. 

As discussed in Sec. IVB 21 the shape of the cell be- 
comes important at very high pressures, and our use of 
the spherical approximation meant that in some cases 
we could only obtain lower limits on the crust thickness. 
Studying more realistic shapes is straightforward in prin- 
ciple, but would require a more demanding multidimen- 
sional numerical solution of the Poisson equation. 

Shell-model corrections can be of order one MeV per 
quark for strangelets of size i? < 5 fm [l^, [2l| , which is 
not negligible relative to our typical enthalpy per quark 
(Fig. [T|) , and may therefore affect our results for the outer 
part of the crust, where we predict strangelets as small 
as 3 fm (Fig. [3D. 

Treating Swiss-cheese phases would require us to in- 
clude curvature energy as well as surface tension. This 
highlights the fact that we treated the interface between 
quark matter and vacuum in a very simplified way, as 
a zero-width interface with a surface tension. However, 
since the quark confinement distance is about 1 fm, the 
interface might well have structure on this distance scale. 
Like the shell-model effects described above, this could be 
relevant to the low-pressure regime, where the strangelets 
can be as small as a few fm. There are even indica- 
tions that when such physics is taken into account, the 
CFL phase may undergo some degree of charge separa- 
tion [2^, [2^ , raising the possibility that there might be 
some sort of crystalline crust on quark stars made of CFL 
quark matter. 
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